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A microscopic analysis of the non-dissipative force F„d acting on a line vortex in a type-II 
superconductor at T = is given. All work presented assumes a charged BCS superconductor. We 
. . . first examine the Berry phase induced in the BCS superconducting ground state by movement of the 

' vortex and show how this phase enters into the hydrodynamic action Shyd of the superconducting 

condensate. Appropriate variation of Shyd gives F^^ and variation of the Berry phase term is seen 
to contribute the Magnus or lift force of classical hydrodynamics to F„d. This analysis, based on the 
BCS ground state of a charged superconductor, confirms in detail the arguments of Ao and Thouless 
O , within the context of the BCS model. Our Berry phase, in the limit e ^ 0, is seen to reproduce 

• the Berry phase determined by these authors for a neutral superfluid. We also provide a second, 

' independent, determination of Fnd through a microscopic derivation of the continuity equation for 

the condensate linear momentum. This equation yields the acceleration equation for the superflow 
and shows that the vortex acts as a sink for the condensate linear momentum. The rate at which 
momentum is lost to the vortex determines Fnd in this second approach and the result obtained 
^ . agrees identically with the previous Berry phase calculation. The Magnus force contribution to F^d 

1/^ ' is seen in both calculations to be a consequence of the vortex topology and motion. 

O 74.20.-z, 03.40.Gc, 74.60.Ge 

(N 

Already in the phenomenological/macroscopic, models of vortex dynamics in type-II superconductors due to 
, Bardeen-Stephen (BS) and Nozieres-Vinen (NV)El, the form of the non-dissipative force Fnd acting on the vortex 
' is controversial. This force is the result of the vortex's interaction with an applied magnetic field Hext, an electric 
field E due to the vortex motion, and the surrounding condensate of superconducting electrons. The disagreement 
centers on whether the vortex feels the lift or Magnus force of classical hydrodynamics as a consequence of its motion 
C through the superconducting condensate. In the BS model, the non-dissipative force is due strictly to the Lorente 
I force pshuj{-Vs x z)/2; while in the NV model, the Lorentz force is supplemented by the Magnus force —psmlCvs x za. 
(~| ' In a very interesting paper, Ao and Thoulesso have returned to this controversy arguing that the correct form for 
O ■ is the NV-form, and that the Magnus force contribution to it is a manifestation of a Berry phase induced in the 

many-body ground state due to the vortex motion. They provide a calculation for a neutral superfluid and argue 
^ that the same scenario will also apply for a charged superconductor. Given that the BCS model of superconductivity 
provides a highly successful microscopic description of the dynamics of a charged superconductor, it would be very 
. _ interesting to see if F^^ can be determined using this model of a charged superconductor (together with the starting 
■ assumptions common to BS and NV, see below) . In this Brief Report we report the results of two such calculations. 
" " ' A detailed presentation and discussion of these calculations will be reported elsewheraj. In the first calculation we 
determine F„£; by working with the BCS superconducting ground state in the case where a vortex is present. This 
state is first constructed and the Berry phase induced in it by the vortex motion is determined. We then show how 
this Berry phase enters into the action describing the hydrodynamic degrees of freedom of the superconducting con- 
densate. Variation of this action with respect to the vortex trajectory gives Fnd and the result found is seen to take 
the NV-form. In the second calculation we give a microscopic derivation of the acceleration equation for the superfiow. 
Together with the expected contributions off the vortex due to spatial variation of the chemical potential, and the 
electric and magnetic fields present, we also find a singular term arising from the vortex topology which describes the 
disappearance of linear momentum into the vortex. The rate at which this momentum is disappearing gives Fnd and 
is found to agree identically with the result of the Berry phase calculation. We stress that the two calculations are 
independent of each other, and each shows that the Magnus force contribution to Fnd arises as a consequence of the 
vortex topology. 

We make use of the Bogoliubov equation to treat the superconducting dynamics. The gap function takes the form 
A(r) = Ao(r) exp[— 10] in the presence of a line vortex with winding number uj = —1 (in cylindrical coordinates (r, 
9, z) centered on the vortex). As in the models of BS and NV, we: (i) assume T = 0; (ii) will approximate the non- 
local character of BCS superconductivity by a local dynamics; (iii) assume H^^ < Hext so that vortex-vortex 
interactions can be ignored and attention can focus on a single vortex; (iv) assume a clean type-II superconductor 
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so that pinning effects can be ignored; and, (v) set h = m = c = L-unless otherwise stated. The sohitions of the 
BogoHubov equation in the presence of a Une vortex are well-knowO and can have positive and negative energies 
relative to the Fermi energy. The superconducting ground state is constructed by occuppying the negative energy 
states. The charge conjugation degree of freedom for the two-component Nambu quasi-particle (NQP) is labeled by 
2s z, and the operator that creates a negative energy NQP is 7„| (where n labels the energy spectrum). Thus, the 
ground state in the presence of a vortex is 

\BCS)^Y[jr.i\0) . (1) 

n 

7„lj-jdepends linearly on the (complex conjugate) of the comuonents of the solutions of the Bogoliubov equation (u„, 
«„)□. Adiabatic motion of the vortex generates a Berry phaseQ 0„ in the solutions {un, ««). Consequently, jni inherits 
the phase — </>„ which, from eqn. (|l|), causes the ground state to develop the Berry phase T ^ — J^n 'f'n- Because the 
electrons are electrically charged, one must use the gauge-invariant form of the Berry phaseQ 



f dT{En\t^ + ^Ao{T)\E,,) 

Jo dr n 



(j)n{t) is calculated using the solutions of Rcf. |^, from which one can then obtain the ground state Berry phase F. One 
finds 

F = y drd^x Qro • Vr„0 - ^A^ , (2) 

where Tq is the vortex trajectory, and we work per unit length of the vortex. We see that our result reproduces the 
Berry phase obtained in Ref. ^ for a neutral superfluid in the limit where e — > 0. 

We now show how the ground state Berry phase F enters into the action describing the hydrodynamic degrees of 
freedom of the condensate. We begin with the vacuum-to-vacuum transition amplitude for the system of electrons 
which can be written as a path integral quadratic in the fcrmion fields via a Hubbard-Stratonovitch transformation 



J V[A]V[A*]{vac;A{t = T)\Us{T,0)\vac;A{0)) . 



Here Ua{T, 0) = T(exp[— i jj' drHeff]); Heff — Hf + L^m + Lc\ Hf is the usual BCS Hamiltonian in the presence 
of a 4-potential {Aq, A); Lem is the Lagrangian for the induced electric and magnetic fields (E, H — He^:*); and Lc is 
the condensation Lagrangian with density \A\'^ /2g. The action for the condensate S = Sq + Shyd is given by 

e-*('^''+^'''"') = (mc; A(r)|t/A(r,0)|wac; A(0)) . (3) 

^0 is the action for the bulk degrees of freedom of the condensate; Shyd is the action for the hydrodynamic degrees 
of freedom; and terms in S containing derivatives of the gap function higher than second order are suppressed. By 
factoring Ua{T, 0) in eqn. (|^) into a sequence of infinitesimal propogations, and appropriately inserting complete sets 
of instantaneous energy eigenkets {\En{tk) )}, evaluation of the matrix element in eqn. (|^) boils down to consideration 
of propogation over an infinitesimal time interval. Spatial translational invariance, which follows from the assumed 
absence of pinning sites, insures that \vac; A(0) ) evolves into the instantaneous ground state \BCS{t) ) of Heff{t), so 
that the relevant matrix element is ( BCS{t + e)|[/A(t)(i -|- e, t)\BCS(t) ). One findsQ 

{BCSit + e)\U/:^(t){t + e,t)\BCS{t)) = e'^'{BCS{t)\e-'"'"^*'^'\BCS{t)) , (4) 

where F is the Berry phase developed in \BC'S{t) ) due to the vortex motion. The remaining matrix element on the 
RHS of eqn. (Q) can be evaluatedQ; and the contribution from all infinitesimal time intervals summed. This yields the 
following result for the hydrodynamic action 
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"^P'-^l + iV(0)io' + ^ I (H - - 



in which the ground state Berry phase F appears as a consequence of the adiabatic motion of the vortex. Here 
Vs — — (?i/2m)[V0 + (2eA)/(?ic)]; is the gap phase; A^(0) is the electron density of states at the Fermi level; 
^0 = gAq + {h/2)dt4>] and Ti, m and c have been re-instated. Appropriate to the scenario of an external current 
passing through a thin superconducting film in the flux-flow regime, we assume the superflow is a combination of an 
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applied superflow v = (h./2m)\/f3 and one that circulates about the moving vortex with velocity Vdrc = —{'h/2m)V0. 
The terms in Shyd linear in Vro6' describe the coupling of the vortex to the applied superflow v; to the electric and 
magnetic fields via (Aq, A); and to the superconducting electrons via the Berry phase T. Variation of the coupling 
terms with respect to ro gives the non-dissipative force 

F„, = ^(v-ro)xz + 0(CoVA2) , 

where is the zero temperature coherence length, and A is the London penetration depth. Our result for 'Pnd is 
identical to the result found by Ao and Thoulesa3 in the case of a neutral superfluid, and which they argued would also 
be true for a charged superconductor. In this first calculation we have considered the case of a charged superconductor 
explicitly (within the context of BCS superconductivity) and found that the Berry phase generated in the SCS ground 
state is responsible for producing the Magnus force contribution to F„d as argued by Ao and ThoulessH, and that 
F„d is given by the NV-result. We go on now to the second independent calculation of ¥nd- 

Our starting point (again) is the Bogoliubov equation for the case where a line vortex with winding number 
Lo = —1 is present. We transform the Bogoliubov Hamiltonian using the unitary operator U — exp[i0(T3/2] to obtain 
Hboq = o'si (*V — (T3Vs)^/(2) — Ef] + AofTi. Here the {ai} are the 2x2 Pwli matrices; Ef is the Fermi energy; and 

= — (1/2)V6' — eA {h — m = c = 1). We make an eikonal approximations for the Bogoliubov equation eigenstates 
(p — exp[iq ■ r](j)' , where |q| = kf and (/)' varies on a length scale L ^ kj^ . To first order in gradients, this gives 
Hsog = cr3[ — q ■ (iV — (T3Vs)] + Aqcti, from which we obtain the gauge-invariant second quantized Lagrangian 



£(q) = *t 



^dt + (73 ( ^dt - eAo ] + cr3q • (iV - (73^^) - Aqcti 



We see that the eikonal approximation made for the eigenstates of Hsog near the Fermi surface in terms of wavepackets 
with mean momentum Pf(\ has led to the separation of the 3+1 dimensional NQP dynamics into a collection of 
independent 1 + 1 dimensional subsystems labeled by directions along the Fermi surface q and which we will refer 
to as q-channels. By construction, both positive and negative energy eigenstates (viz. above and below the Fermi 
surface) carry a mean momentum p/q. Positive energy quasiparticles in this channel carry (mean) momentum Pf(\ 
(right-goers, V"!?)! while positive energy quasiholes have (mean) momentum — p/q (left go-ers, and spin indices 
have been suppressed. The adjoint of the NQP field operator in this channel is ^'^(x) = ('0]j(x;q) VL(x;q)). The 
Noether current associated with the global phase transformation ^I^q exp[— ix]^'q can be written in a psuedo- 
relativistic notation as = ^7'^\E'. Here /i = 0, 1; x° = t, a:"'^ = q • x; 7° = ui, 7^ = —ia2\ and ^ = ^[''''7°. One 
can then write the density of linear momentum (in the q-channel) as gi(a;; q) = Pfiiij'^{x)] and the associated stress 
tensor as Ty (x;q) = Pf^:i^.ji^{x)■ Taking the expectation value of these operators with respect to the q-channel 
ground state |?;ac)q, and summing over all q-channels gives the ground state density of linear momentum gi(x) in 
the condensate and its associated stress tensor Tij{x) 

g,{x) = kjY^ f ■^ci,{vac\f{x;q)\vac)q ; T,j{x) = kjY^ /" ^ q^q^ ( i;ac|j^(a;; q)|i;ac )q ; (5) 

where a is the spin index (±). The continuity equation for the condensate linear momentum is then 

dtgt + djT.j ^ kjY^ /" ^ q, ('uac|9^j^|t)ac)q . (6) 

The matrix element appearing in cqn. ^ docs not vanish, signaling that the condensate linear momentum is not 
conserved (not surprising since the condensate is not isolated). Wc will see shortly that, together with the expected 
source terms due to gradients in the chemical potential, and from the electric and magnetic fields; there will also be 
a source term whose origin lies in the vortex topology and which keeps track of the rate at which linear momentum 
is disappearing into the vortex. This topological term will thus give F„d in this second approach. Details of the 
calculation of M = (vac\df_ij'^\vac)q are given in Ref. |^. The result is M = {e^'^ F^^^) / in , where e"^ = 1; i^^^ = 
di_iA^ — d^Afj,; and Aq — cAq — (1/2)9(6'; Ai = q • v^. Inserting this result for M into eqn. (^ gives 

^tg^ + ^JT,J^Co^ -^[do,d,]0 + eEA . (7) 
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Here Cq = fc^/37r^ is the particle density in the normal phase for the case where the chemical potential equals the 
Fermi energy; and h has been re-stored. The first term on the RHS of eqn. (^ is non-vanishing due to the non- 
trivial vortex topology. The local expression of this topology, appropriate for a vortex with winding number uj, is 
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[dxTdy\(j) = 2ttui5{x — xo)5{y — yo), where is the gap phase and rp = {xn, yo) is the position of the vortex. The 
LHS of eqn. can also be evaluated using eqn. ^ along with the resullH {vac\j'^\vac)q = Anj^-K. The results 
are ~ Co (vs)^; ~ Co {^dtO — eAo) Sij. Making use of these results in eqn. (|^); together with the Josephson 
equation (hdt4>)/2 = —fioj where /Xq is the chemical potential in the vortex rest frame which can be written as 
/io = M + + eAo (/i is the chemical potential in the lattice frame and m — 1) gives finally 

^^ = -VM + eE + ev, xB-^(v, -ro)xz(52(r-ro) . (8) 

We see that the continuity equation for the condensate linear momentum has yielded the acceleration equation for the 
superflow. We find the expected source terms related to the hydrodynamic pressure (VP = p^V/i), and the electric 
and magnetic fields. We also see that linear momentum is disappearing from the condensate into the vortex at ro(t) 
at the rate {pshuj/2){vs — r'o) x z per unit length so that 

Fnd = - ro) X z , 

in agreement with the Berry phase calculation. Our result is also consistent with the calculation of NV in Rcf. 
These authors showed that the first 3 terms in eqn. (|^) lead to a flux of linear momentum in towards the vortex at 
a rate (pshuj)/ (2)(ys — r'o) x z which is exactly the rate at which we find it appearing on the vortex, indicating that 
linear momentum is conserved in the combined condensate-vortex system. 

In this paper we have provided two independent microscopic calculations of the non-dissipative force F„d acting 
on a line vortex in a type-II superconductor at T = 0. Both calculations yield the NV-form for this force Fnd = 
{pshuj /2)(y — r'o) x z. The first calculation (inspired by earlier work of Ao and Thouless which determined F„rf via 
a Berry phase analysis appropriate for a neutral superfluid, and which they argued would also be true for a charged 
superconductor) shows that the arguments of Ao and Thouless are fully borne-out in the context of the BCS model for 
a charged superconductor. The second calculation (which does not rely on Berry phases) examines the flow of linear 
momentum in the condensate. The continuity equation for this linear momentum is shown to: (i) yield the acceleration 
equation for the superflow; and (ii) to contain a sink term indicating the disappearance of linear momentum into the 
vortex. F„d follows in this second approach from the rate of momentum loss to the vortex. The result obtained is 
the NV result and the Magnus force contribution to F„d is seen to be a consequence of the vortex topology. 

Note Added: Two preprints have appeared since this work was completed (M. Stone; Aitchison et. al.ll3) which also 
find a gauge-invariant contribution to the hydrodynamic action first order in time derivatives of the gap phase. 

I would like to thank Ping Ao for interesting me in this problem and for helpful discussions; Michael Stone for 
interesting comments and discussions, particularly with regard to the issue of gauge invariance; T. Howell III for 
constant support; and NSERC of Canada for financial support. 
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